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ABSTRACT 

We combine predictions for several hierarchical cosmogonies with observational evidence on damped 
Lya systems to establish a correspondence between the high redshift galaxy population and the properties 
of damped Lya systems (DLAS). We assume that high redshift galaxies and damped Lya systems are 
hosted by the same dark matter halos and require consistency between the predicted halo space density, 
the rate of incidence and the velocity width distribution of damped Lya systems, and the observed galaxy 
luminosity function at the bright end. We arrive at the following results: (1) predicted impact parameters 
between the damped absorption system and the luminous parts of the absorbing galaxy are expected to 
be very small (0.3 - 1 arcsec) for most galaxies; (2) luminosities of galaxies causing damped absorption 
are generally fainter than = 25 and damped Lya systems are predicted to sample preferentially the 
outer regions of galaxies at the faint end of the galaxy luminosity function at high redshift. Therefore, 
DLAS should currently provide the best probe of the progenitors of normal present-day galaxies. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — quasars: absorption lines 



1. INTRODUCTION 

The physical conditions inferred from the absorption 
features caused by high redshift damped Lya absorption 
systems (DLAS) are, in several aspects, similar to those 
in the interstellar medium of present-day galaxies. It has 
therefore been suggested to identify DLAS with the pro- 
genitors of such galaxies (e.g. Wolfe 1988). At low redshift 
DLAS show a wide variety of morphologies (Le Brun et al. 

1997) . At high redshift, however, few have been detected 
in emission (M0ller & Warren 1995; Warren & M0ller 1995, 
Djorgovski et al. 1996, Djorgovski 1997, M0ller & Warren 

1998) and the nature of galaxies causing the absorption 
remains unclear for the majority of DLAS at high red- 
shift. The main source of information on the nature of 
DLAS comes from the associated metal absorption which 
probes the kinematics and the chemical enrichment his- 
tory of the mainly neutral gas in these systems. Motivated 
by the characteristic asymmetric shape of the absorption 
profiles of low ionization species Wolfe and collaborators 
(e.g. Wolfe 1986) have suggested that high-redshift DLAS 
are large rapidly rotating discs akin to present-day spiral 
galaxies (but see also Ledoux et al. 1998). In this pic- 
ture the generally low observed metallicity of high red- 
shift DLAS (Pettini 1994, Lu 1996, Pettini 1999) is taken 
as evidence that these are still chemically young (Lindner, 
Fritze-von Alvensleben & Fricke 1999). 

The conjecture of DLAS being large rotating discs is, 
however, at odds with the velocity width distribution, the 
size and the total cross section of rapidly rotating discs pre- 
dicted by hierarchical cosmogonies (Haehnelt, Steinmetz 
and Ranch 1998, hereafter HSR98). Hydrodynamical sim- 



ulations (Katz et al. 1996; Haehnelt, Steinmetz & Ranch 
1996; Gardner et al. 1997ab; Ranch, Haehnelt & Steinmetz 
1997) show that gas condensations leading to damped Lya 
absorption do indeed occur in such scenarios. However, the 
hierarchical cosmogonies predict that a significant fraction 
of the cross section for damped absorption is contributed 
by gas in halos with circular velocities as small as 50 km s^^ 
(Ma & Bertschinger 1994; Mo & Miralda-Escude 1994; 
Klypin et al. 1995, Kauffmann 1996, Mo, Mao & White 
1999; Gardener et al. 1999), where we define the circular 
velocity as Vc = \/ GMjr at a radius where the overden- 
sity is equal to 200. In HSR98 we have shown that the 
shape of the absorption profiles is not unique but can be 
equally well produced by merging proto-galactic clumps 
(cf. Nulsen, Barcon & Fabian 1998; McDonald & Miralda- 
Escude 1999, Mailer et al. 1999). In HSR98 we further 
showed that the observed velocity width distribution of 
the absorption features of low ionization species can be 
reproduced in a standard cold dark matter (SCDM) cos- 
mology. At the same time hierarchical cosmogonies such 
as SCDM have been shown to reproduce the basic prop- 
erties of the observed population of high-redshift galaxies 
(Adelberger et al. 1998; Steidel et al. 1998; Baugh et al. 
1998; Bagla 1998; Jing & Suto 1998; Haehneh Natarajan 
& Rees 1998; Jing & Suto 1998; Contardo, Steinmetz & 
Fritze-von Alvensleben 1998; Steinmetz 1998; Kauffmann 
et al. 1999; Mo, Mao & White 1999). Here we estab- 
lish a link between emission and absorption properties of 
high-redshift galaxies. In section 2 we briefly test if our 
previous findings depend on cosmology. In section 3 we 
predict the luminosity and impact parameter distribution 
of high-redshift galaxies responsible for DLAS. Section 4 
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Table 1 

Model parameters: as is the rms linear ovcrdcnsity in spheres of radius 8h^^ Mpc and F is a shape parameter for CDM-like spectra, h is the 
Hubble constant (in lOOkms"^) and Clo S'lid CIa a-i'e the total energy density and that due to a cosmological constant. For the other 

quantities see text. 
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contains our conclusions. 

2. THE VELOCITY WIDTH DISTRIBUTION OF DLAS 

In hierarchical CDM-like cosmogonies DLAS arise nat- 
urally from the cool gas that accumulates at the center 

of dark matter (DM) halos (e.g. Katz ct al. 1996, Gard- 
ner et al. 1997a/b). Absorption features of low ionization 
species like Sill are generally believed to be good tracers 
of the motion of this gas in a gravitational potential well 
(Prochaska & Wolfe 1997, 1998). One major characteristic 
of these absorption features is their overall width which in 
HSR98 we have shown to be correlated with the circular 
velocity of the DM halo hosting the DLAS. However, for a 
halo of given circular vc^locity a statistical distribution of 
velocity widths p((wwid| '^'c) arises due to different orienta- 
tions of the linc-of-sight and different dynamical statc!s of 
the gas of the DM halos. Halos with small circular veloci- 
ties are believed to lose most of their gas due to feedback 
effects. The velocity width distribution of DLAS as probed 
by low ionization species can thus be written as 

/•OO 

Pdamp(Wwid) = ^ P{Vwid\Vc) X Pdamp(t^c) dv^. (1) 

In HSR98 we used simulated absorption profiles for a set 
of haloes with 50kms~^ < Vc < 250kms~^ to investigate 
p(^^widl^^c) for a SCDM model and found that the distribu- 
tion does not depend explicitly on Vc and depends little on 
redshift. We therefore assume p(wwidbc) to be a function 
of Wwid/i'c only. In order to check if there is a dependence 
on cosmology we repeated our analysis for a ACDM model 
(see table 1 for the assumed model parameters and HSR98 
for details of the simulations). Figure 1 shows p{vv,id/vc) 
averaged over different halos and lines-of-sight. There is 
little difference between the two cosmological models. 

Once the cross-section weighting for damped absorp- 
tion is known the velocity width distribution is read- 
ily calculated from the distribution of circular velocities 
i^damp (vc) oc fJdamp X Pps(^'c)- K is Currently difficult to 
infer this cross-section weighting reliably from numerical 
simulations. Gardner et al. (1997) e.g. found that the 
cross section of DM halos scales as a power-law of the cir- 
cular velocity as v^'^^^. In their later work they found 



a much shallower scaling v^a-li (Gardner et al. 1999). 
This cross-section weighting is likely to be sensitive to the 
energy and momentum input due to supernovae which is 
difficult to include properly in these simulations. 
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Fig. 1 . — The distribution of velocity widths of absorption features 
of low ionization species as a function of v^\^/vc for two different 
cosmogonies. Model parameters are given in table 1. 

Here we take a different approach and use the observed 
velocity width distribution to determine the mean cross 
section weighting for damped absorption. We also choose 
a power-law parameterization 

<7damp(t^c) = 7r(rl^p)^(t;c/200kms-i)'5. (2) 

Using the Press-Schechter formalism we calculate Pps(wc) 
and obtain the velocity width distributions shown in Fig- 
ure 2. The observed distributions are well fit with /? ~ 
2.5 — 3 (see table 1) and iJniin = 50kms^^. Note that 
there is little evolution with redshift in the observed as 
well as in the predicted distribution. The values of beta 
are close to the old value but quite different to the new 
value of Gardner et al. (1997,1999). 

3. PREDICTING EMISSION PROPERTINGS 

3.1. The luminosity distribution of DLAS 

Recently, it has become possible to determine the lumi- 
nosity function for the population of star-forming galaxies 



Haehnelt, Steinmetz & Rauch 



Vol. 999 



15 



I 10 
> 



5 - 



n — I — rrj — I — i — I — I — I — I — I — I — 1 — r 

SCDM 



"T — r~3 — I — I — I — I — I — I — I — I — I — 1" 

ACDM 



T — I — rrj — I — I — I — I — I — I — I — 1— 

OCDM 




~? — I — I — I — I — I — I — I — I — I — r 

tCDM 



too 200 
v,y/kms-' 



100 200 



100 200 
v,.ykms-' 



100 200 
v,ykms-' 



Fig. 2. — The velocity width distribution of the associated absorption of low ionization species in DLAS. The curves are model predictions 
for four different cosmogonies as calculated from equation (1) (parameters see table 1). The crosses are data from Prochaska & Wolfe (dotted: 
high z, solid: low z). 



detected at high redshift (Steidel et al. 1996, Steidel et al. 
1999). The himinosity function at z=3 can be reproduced 
if a simple linear scaling of the luminosity in the TZ band 
(Steidel et al. 1996) with mass of the DM halo is assumed, 

TO^ = - 7.5 X log[i;c/200kms"^] (3) 

(Haehnelt, Natarayan & Rees 1998). The required val- 
ues of are also shown in table 1. The same simple 
model can also reproduce the clustering strength of high- 
redshift galaxies and its decrease with decreasing limiting 
UV luminosity of the galaxy sample (but see Somerville, 
Primack & Faber (1999) for a somewhat different model). 
Note that the linear scaling of the UV luminosity with 
mass does not necessarily imply that the star formation 
rate scales linearly with the mass of the DM halo as the 
dust extinction of the UV luminosity is probably luminos- 
ity dependent (Steidel et al. 1999). If we make the plau- 
sible assumption that high-redshift galaxies are hosted by 
the same population of dark matter halos that are causing 
DLAS and assume that the linear scaling of the UV lu- 
minosity with mass can be extrapolated to smaller masses 
we can predict the luminosity distribution of DLAS. All 
we need to know is the scaling of the cross-section of DM 
halos for damped absorption with circular velocity which 
we already inferred in the last section from the velocity 
width distribution. The result is shown in Figure 3. 80 
to 90 percent of DLAS are predicted to be fainter than 
= 25.5 (the spectroscopic limit), rather independent 
of cosmology. Some of these may be faint Lya emitters 
(Fynbo, Thomsen & M0ller 1999) but most of them will 
probably be very difficult to identify reliably. 

3.2. The impact parameter distribution of DLAS 

Similarly we can also predict the impact parameter 
distribution (Fig. 3b). The overall rate of incidence, 
dA/"/ dz = 0.2 (Storrie-Lombardi et al. 1996), determines 
the normalization of the cross-section weighting. The cor- 
responding values of r'^^^^ = rdamp('t'c = 200kms~^) are 
given in table 1. The predicted impact parameters range 
between 0.3 and 1 arcsec, rather independent of cosmol- 
ogy. The parameter r^oo = '"20o('i^c = 200kms~^), i.e. the 
radius at which the overdensity is equal to 200 x pcrit is 
also listed in table 1. 



The typical radii of the region causing damped absorp- 
tion are generally about 25 to 50 percent of r2oo of the 
corresponding DM halo. This is about a factor ten larger 
than the expected scale length of a centrifugally supported 
disc if the angular momentum of the gas is due to tidal 
torquing during the collapse of the DM halo (Mo, Mao & 
White 1998). Note that with the assumed scaling of the 
cross section rjamp increases more steeply with circular 

velocity (oc Vc^^) than r2oo (oc Wc). 

4. CONCLUSIONS 

In hierarchical cosmogonies DLAS are mainly caused by 
merging protogalactic clumps hosted by collapsed DM ha- 
los. The observed velocity width of the absorption profiles 
of associated low ionization species is well reproduced for a 
a range of CDM variants which reproduce the present-day 
space density of galaxy clusters. 

By combining the cross-section weighting inferred from 
the observed velocity width distribution with the mass lu- 
minosity relation inferred from the luminosity function and 
clustering properties of high redshift galaxies it is possi- 
ble to link absorption and emission properties of DLAS. 
About 10 to 20% of DLAS are predicted to be brighter 
than = 25.5. Expected impact parameter typically 
range between 0.3 and 1 arcsec but there is a pronounced 
tail of larger impact parameters. The predictions seem 
consistent with the luminosities and impact parameters of 
the small number of DLAS with spectroscopically iden- 
tified emission (see e.g. M0ller & Warren 1998 for an 
overview). In our model these would be drawn from the 
tail of the distribution at the bright end and at large im- 
pact parameter, respectively. The rather faint flux levels 
and small impact parameter predicted for the majority of 
DLAS explains why searches for the emission from DLAS 
have been notoriously difficult. 

The predicted impact parameters are nevertheless about 
a factor three to flve larger than the typical radii of the 
luminous regions of Lyman-break galaxies. Thus, DLAS 
at high redshift should preferentially sample the outer re- 
gions of galaxies at the faint end of the luminosity func- 
tion. This is consistent with the low observed metallicities 
in DLAS and the outer regions of galaxies. In our model 
DLAS should not be expected to show the high metal- 
licities estimated for Lyman break galaxies (Pettini et al. 
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Fig. 3. — Left: The predicted luminosity distribution for the DLAS 
parameter distribution. 

1999); the latter are derived from emission spectra, which 
are dominated by the light from the more highly enriched 
central regions. 

Hierarchical cosmogonies predict Lyman-break galaxies 
with < 25 to become part of bright galaxies in galaxy 
clusters at the present day. The fainter galaxies responsi- 
ble for DLAS should, however, be the building blocks of 
more typical present-day galaxies consistent with the 
findings from chemical evolution models for DLAS. Most 
high-redshift DLAS are predicted to be too faint to obtain 
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spectra even with 10m telescopes unless they are strongly 
magnified by gravitational lensing. The analysis of their 
absorption properties is currently the prime method for 
studying the progenitors of normal present-day galaxies. 
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